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■^ Abstract 

I— I We used the internal conversion {Eq transition) of germanium-72 to indirectly mea- 

Q sure the low energy nuclear recoils of germanium. Together with a reliable Monte 

Carlo package, in which we implement the internal conversion process, the data 

was compared to the Lindhard (/c=0.159) and Barker-Mei models. A shape analysis 

indicates that both models agree well with data in the region of interest within 

c/) 4%. The most probable value (MPV) of the nuclear recoils obtained from the shape 

• ^H analysis is 17.5 keV with an average path-length of 0.014 /um. 
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Understanding detector response to low energy nuclear recoils is imperative to 
the interpretation of experimental results from a detector designed to search 
for WIMPs (weakly interacting massive particles), a dark matter candidate. 



;h Direct detection of low mass WIMPs occurs in the low energy region of detec- 



tors with a threshold down to sub keV. Since the threshold energy represents 
the visible energy in the detector, understanding the ionization efficiency of 
the detector response to low energy nuclear recoils is crucial to calculating the 
recoil energy. An example of the need for this requirement is the claim of exper- 
imental evidence for dark matter by CoGeNT [Ij that has been unverified by 
CDMS II P]. Both experiments use germanium as the target material. Thus, 
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similar results are expected if the detection thresholds for both experiments 
were determined using a standardized ionization efficiency, which accurately 
accounts for all processes that occur at a low energy range. This standardized 
ionization efficiency must also be validated with measurements in combination 
with reliable Monte Carlo simulations. 

Two different approaches can be used for modeling ionization efficiency in 
germanium detectors J3|H] . One model, traditionally used for a number of 
different detector materials, is proposed by Lindhard et. al. [3]: 



l + k-g{ey 



where, Z is the atomic number, A is the mass number, E^ is the recoil energy, 
k = 0.133^2/3^-1/2^ g(e) = 3e°-i^ + 0.7e°-'5 + e, and e = ILSE^Z'^/^. However, 
this model has not been proved accurate at low energies as the theoretical 
derivation has uncertainties in this region [5]. 

Another model designed for low energy interactions in germanium was pro- 
posed by Barker and Mei [1]. This model takes into account the fraction of 
nuclear stopping power that contributes to the ionization efficiency at low 
energies [4j. The Barker-Mei model can be expressed as: 

0.14476 ■^0-697747 

[2) 



' -1.8728 + exp[E0-2ii349]- 

This model is valid for recoil energies, Er, from 1 keV to 100 keV. However, the 
Barker-Mei model has not been experimentally verified; this paper addresses 
that issue. 

A comparison between the two models and existing data was performed [1], 
and Fig. [T] shows both models together with available experimental data. As 
shown in Fig. [l} the two models agree in the low energy region but disagree 
with the available data points in the same region. Further validation of the 
Lindhard and Barker-Mei models, using more measurements and accurate 
Monte Carlo simulations, is necessary. 

Taking exact nuclear recoil measurements can be challenging at lower energies 
when systematic errors are introduced by a variety of sources. For example, 
thermal neutrons and elastic/inelastic scattering will have their own uncer- 
tainties. Utilizing thermal neutrons requires the full absorption of out-going 
gamma rays measured by another detector in coincidence. Without this addi- 
tional measurement, Compton scattering from out-going gamma rays within 
the germanium detector can contaminate the visible energy. In a measure- 
ment of neutron elastic scattering, a Monte Carlo simulation must be incorpo- 
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Fig. 1. The Barker-Mei and Lindhard models of ionization efficiency with the ex- 
perimental data !5l()l7l8l9ll()|ll|12|i:Sj . 

rated to exclude multiple scatters. It is also necessary to precisely measure the 
scattering angle and time-of-flight of the out-going neutrons. With inelastic 
scattering, the Compton scattering due to de-excitation of gamma rays in the 
detector can contaminate the signal. All of these techniques need to be imple- 
mented with an accurate Monte Carlo that reduces systematic errors in order 
to obtain a reliable ionization efficiency. However, current popular simulation 
tools often need to be tuned to simulate inelastic scattering processes. 

Therefore, we desired a simple method to accurately measure nuclear recoils 
in germanium. The Eq transition of germanium- 72 (^^Ge(n,n e)), which is the 
internal conversion process for this nucleus, was chosen. The Eq transition of 
^^Ge(n,n e) is induced when neutrons inelastically scatter off a ^^Ge nucleus. 
After this collision, the ^^Ge nucleus is left in an excited state. When the nu- 
cleus returns to ground state, it does not produce a gamma-ray (as is common 
for other nuclei), but interacts electromagnetically with the inner shell elec- 
trons and causes one to be emitted from the atom [13]. The physics process 
for the internal conversion of ^^Ge can be generalized as: 



n+"Ge 



n' +^2 Ge* -^^2 Q^ + g- 



+ X — ray, 



(3) 



where ^^Ge* denotes the excited O"*" state. The total detectable energy from 
the '''^Ge*(0"'") to '^^Ge(0"'") transition is 691.6 keV, which includes energy from 
the X-ray. This energy is well within the visible range of traditional detectors, 
and both electron and X-ray can be detected with a consistency of approxi- 
mately 100%. The characteristic energy 691.6 keV is superimposed with the 



nuclear recoil of a ^^Ge nucleus to form a quasi-triangular shape, which is dis- 
tinguishable from other processes in the spectrum. Thus, we can observe low 
energy nuclear recoils without taking low energy nuclear recoil measurements. 
The quasi-triangular shape, starting at 691.6 keV, has been studied in great 
detail p3|20f21f22|l23] . The tail of the quasi-triangular shape is a combination 
of the 691.6 keV energy deposition and nuclear recoil energy due to neutron 
scattering. 

Depending on the incident neutron energy and scattering angle [15], nuclear 
recoil energy can be expressed as: 

where E^ is the nuclear recoil energy, E„ is the neutron energy, M is the nucleus 
mass, m is the neutron mass, and 9 is the scattering angle between the incident 
neutron and the recoil nucleus. The quasi-triangular shape is created by adding 
the energy from the nuclear recoils to 691.6 keV energy deposition as: 

E = mi.QkeV + fEr, (5) 

where E is the observed energy and e is the ionization efficiency. The nuclear 
recoil energies can be determined using a Monte Carlo simulation with applied 
ionization efficiency, if an agreement with the measurements is obtained. 

To validate the ionization efficiency models proposed by Lindhard and Barker- 
Mei,a Geant4.9.2-based Monte Carlo simulation package, corrected for the in- 
ternal conversion processes [16], was used. This simulation package was verified 
using a well-calibrated ^°Co radioactive source with the same experimental 
setup. 

Utilizing the measurements with a substantiated Monte Carlo simulation, we 
compared the unique quasi-triangular shape induced by the 691.6 keV elec- 
trons and nuclear recoils in the data with Monte Carlo simulations for the 
two models. Shape analysis was used to verify the quasi-triangular shape of 
spectra in the data and Monte Carlo simulations by analyzing data points in 
the region of interest. We found a good agreement between measurements and 
the Lindhard (with k = 0.159) and Barker-Mei models. 

In this paper, we corroborate the ionization efficiency models proposed by 
Lindhard et al and Barker-Mei with a neutron induced Eq transition for ^^Ge. 
The experimental design is discussed in Section 2 followed by data analysis 
in Section 3. The Monte Carlo simulation is demonstrated and validated in 
Section 4, and the comparison with data described in Section 5. Section 6 
summarizes our results. 



2 Experimental Design 



The germanium detector used in our experiment was an old coaxial detector 
from Princeton Gamma Tech, model RGllB/C [17]. Its linearity of energy re- 
sponse between consecutive calibrations was within 0.35%. However, the full 
width at half maximum (FWHM) was 7.1 keV at the 1173 keV ^°Co peak, 
a factor of two worse in energy resolution than a new germanium detector. 
Nevertheless, the detector was sufficient for measuring neutron induced inter- 
nal conversion. We used an ^^^Am-^Be (AmBe) source with neutron energies 
ranging from ~1.0 to 11.2 MeV at a frequency of 100 Hz [18] . 

An AmBe neutron source produces neutrons in four discrete groups, no, n-i, ^2, 
and na, which populate the ground state, the 4.443 MeV level, the 7.65 MeV 
level, and the 9.64 MeV level of the ^^C product nucleus, respectively |2^|l25] . 
The Til group neutrons, accompanied by 4.443 MeV gamma rays, dominate 
the production of neutrons in an AmBe source [2^f25] . This feature allows 
us to set up a coincidence measurement using sodium-iodide (Nal) detectors 
(Bicron model number 3M3/3 [I9]) with a threshold of above 1 MeV. Utilizing 
this coincidence method, the Nal detectors will measure gamma rays while 
neutrons are being detected with the germanium detector. This coincidence 
pattern required that the Nal and Ge detectors both trigger in order for an 
event to be recorded, suppressing random background events generated by 
gamma rays from surrounding materials. We took approximately 22 days of 
data with three Nal detectors and 13 days of data with two Nal detectors. 

Due to the low neutron emission, the source was placed on the center of the 
Ge detector cap and held in place by electrical tape. When using three Nal 
detectors, they were placed on the right, the left, and directly in front of 
the germanium detector. When using two Nal detectors, they were placed 
on the right and on the left of the germanium detector. Fig. [2] shows the 
experimental set-up. The data acquisition was performed using a National 




Fig. 2. Experimental setup for Ge and Nal detectors. Left: Three Nal detectors. 
Right: Two Nal detectors. 



Instruments PXI-1031 system [26] and Igor Pro 4.07 software [27]. Each run 
lasted approximately 4.5 days. The data from each individual run was added 
consecutively off-line using analysis code from the Root software package [28j . 
Evidence of the 691.6 keV Eq transition peak can be seen in Fig. [3j 
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Fig. 3. Data taken with AniBe source after 35 days. 



3 Experimental Results 



3.1 Data Analysis 



After 35 days, the number of events in the region of interest (670 ~ 760 keV) 
was approximately 500. This is consistent with the known neutron emission 
from the AmBe source and coincidence method, which suppresses random 
background events from gamma rays. Because of the small sample size of 
valid data in the region of interest, the bin size in the histogram is set to 6 
keV in order to mitigate statistical fluctuation. As shown in Fig. [3} there are 
several peaks near the region of interest. We identified these peaks in order to 
help understand the processes that took place within our detector as well as 
possible sources of contamination. 



The first peak can be fitted using a standard Gaussian distribution and linear 



background distribution. This is given by: 



Po ■ exp 




P3 + piE, (6) 



where po = 37 is the normahzation constant, pi = 508 keV is the center 
value, p2 = 4.7 keV is the Gaussian width, p^ = -65 is a constant, p^ = 0.2 is 
the slope, and E is the energy in keV. The peak at 511 keV is mainly from 
the annihilation of e"'"e~ pairs induced by cosmic rays passing through the 
surrounding materials. The fitted central value of 508 keV is slightly lower 
than the expected 511 keV, but still within the margin of error for the bin 
size. 

The remaining peaks can be fitted using a Moyal distribution and linear back- 
ground. The Moyal distribution used is [29] : 



* =/ ^^^P [- (^ (^ - ^-P-) + ^^P [^ (^ - ^-P-)])]' (7) 



where i? is a constant and E^pv is the most probable value of energy deposition 
in the detector. R and E^pv are physically significant parameters. 

The peak around 560 keV, is likely caused by events from '^^Ge(n, n'7) inelas- 
tic scattering. Because of the small number of events, the fitting function was 
only partially accurate. The second peak at around 600 keV has the fitted 
parameters, R = 0.06 and Empv = 609 keV. This peak is likely the combi- 
nation of inelastic scattering from ^^Ge(ra,n'7) and the neutron capture on 
^^Ge, ^^Ge(77,, 7). Finally, the peak of interest, 691.6 keV, exhibits the internal 
conversion of ^^Ge(n,n'e). The fitted parameters are, R = 0.06 and E^pv = 
705 keV. All fits in the region of interest are shown in Fig. |4j 

We also took a background spectrum using the coincidence method to identify 
any random coincidence. This spectrum is also shown in Fig. |4j The back- 
ground data was taken using two sodium iodide detectors (see right of Fig. [2]) 
for 4.6 days. 

After the background spectrum was collected, we subtracted it from the AniBe 
data sets as shown by Fig.[5j The 511 keV peak was significantly reduced, veri- 
fying our conjecture on the main origin of this peak as induced pair production 
from cosmic rays in the surrounding materials. The subtracted data set was 
also fitted using a Moyal distribution and linear background. It was found that 
the most probable value of energy deposition in the detector for the Eq tran- 
sition was 696 keV. Thus the percent difference of the AmBe data with and 
without background subtraction is 1.3%. This difference is within the error 
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Fig. 4. Region of interest with all peaks fitted, and the coincidence background 
spectrum. 
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Fig. 5. Region of interest with the background subtracted from the original data 
files. 



margin of a histogram with 6 keV bin size. 



3.2 Interpretation of Results 



Electrons and nuclear recoils that travel through germanium with a high mo- 
mentum lose energy by exciting and ionizing the germanium atoms. The av- 
erage amount of energy lost can be calculated with the Bethe-Bloch equa- 
tion [30]. However, the energy transfer is not a continuous process. It occurs 
through random collisions during which various amounts of energy can be 
transferred. The energy loss in the detector can be described by the Moyal 
distribution, Eq. [7| as shown in Fig. |4} 

R is the reciprocal of the density, p, times the average path-length of charged 
particles in the detector, d, times the parameter, K, given by the Bethe-Block 
formula: 

K ■ p ■ a 

where K = AnNanieC^r^z^ ^^2 is related to the stopping power, with constants 
Na, Avogadro's number; nie, the mass of electrons; c, the speed of light; Ve, 
the Bohr radius; z, the electron charge; Z, atomic number of target; A, the 
atomic mass number of target; and (3, the speed of charged particles divided 
by c. The density of germanium is p = 5.323 g/cm^. 

Since electronic stopping power is different from nuclear stopping power, the 
value of R is very different for pure electronic recoils and nuclear recoils. 
Thus, the value of R can be used to derive the average path-length for a 
given electronic recoil or nuclear recoil. However, the value of R from our 
measurements is a combination of electronic and nuclear recoils, and it cannot 
be used to directly determine the average path length for either. Nevertheless, 
we can use the most probable value of energy deposition, which is related to 
the stopping power multiplied by the average path length, to determine the 
average path length. 

Given Empv = 696 keV, obtained from the fitted function in Fig. |5| is a con- 
volution of the 691.6 keV energy deposition and the nuclear recoils, then we 
contend that 696 keV = j^ -d ■ p, where j^ is mass stopping power in keV 
cm^/g. Since there are both electronic and nuclear recoils, we can rewrite this 
as: 



where |g = 1301 keV cmVg |3l] and jg ■ 4 ■ P = 691.6 keV; thus 4 = 
691.6 /{-^ ■ p) = 0.1 cm, which is the average path length in germanium for 



electrons with an energy of 691.6 keV. Therefore, we have 

dE 
696 fce1/ = 691.6 keV + -—-dn-p (10) 

dXn 



and 



dW 

AA keV = —— -dn-p. (11) 



dXn 



From Eq.(ll), we can conclude that the most probable nuclear recoils that 
we measured in the detector have an average of 4.4 keV electronic equivalent 
energy, which corresponds to approximately 17.5 keV nuclear recoil energy 
from ionization efficiency in the Barker-Mei model. 

To determine dn, we can use j^ = 2341720 keV cm^/g for a nuclear recoil of 
17.5 keV 0]. Thus, d„ = 17.5 keV /(^ ■ p) = 0.014 pm is the average path 
length in germanium for 17.5 keV nuclear recoils. This determines the range 
of low energy nuclear recoils in a germanium detector. 



4 Monte Carlo Simulation of En 



An accurate Monte Carlo simulation is needed to determine the validity of the 
Lindhard and Barker-Mei models when compared to collected data. Two cru- 
cial steps were taken before creating a Monte Carlo simulation that would de- 
termine nuclear recoil energy: 1) the Eo transition in Geant4.9.2 [16j was fixed 
and 2) the Monte Carlo simulation was verified with a well-known gamma-ray 
source. 



4-1 Fixing the Eq transition in Geant4.9.2 



By studying the inelastic scattering processes that contribute to internal con- 
version, we found that the Eo transition code is included in Geant4.9.2, but is 
missing neutron data for ^^Ge. Specifically, Geant4.9.2 does not provide 7/e 
transition data in each energy level (Data Type = 12) and has no cross section 
data corresponding to the Eq transition. 

In order to make the internal conversion process available in our Geant4.9.2 
simulation, we created a 7/e ratio in the database for ^^Ge (see our correction 
in Ref. [22] )• In addition, we calculated all the cross-sections of the ^^Ge(n, n') 
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Fig. 6. The orbital e transition ratio of ^^Ge in its excited states. 

reaction using TALYS [33], a reliable software for the simulation of nuclear re- 
actions. We then converted the transition cross-section data into the transition 
ratio as shown in Fig. [6} 

The independence of this ratio can be cross-checked by looking at the reac- 
tions in different energy scales. Eventually, the TALYS data was converted to 
Geant4 format and added to the bottom of the file, 

G4.NEUTR0NHPD AT A/ Inelastic/ F01/32_72_Germanium, 

where it could be used in the internal conversion process (see the supplemental 
data used in Geant4.9.2 from Ref. f32]). 

During the Eq transition, a characteristic X-ray, or Auger electron, is produced 
simultaneously with the conversion electron. This process was not included in 
the Geant4 code. However, since this is a complicated process that involves 
binding energy, it would require editing for all the elements in the Geant4 
package. For the sake of simplicity, we combined the X-ray and conversion 
electron together, since they both contribute to the total energy deposition 
(in a solid germanium detector only). A more general solution will be required 
when considering other materials. 
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Fig. 7. Simulated geometry for Ge and Nal detectors. Left: Cross section of Ge de- 
tector. Right: Ge detector (Blue), two Nal detectors (Red) and Lead bricks (Gray). 

A simulation was created from the Geant4.9.2 and G4NDL3.13 packages that 
includes the additional transition ratio data. For further information, readers 
can refer to Ref. p7] for the detailed geometry and experimental framework. 
After the missing data was added, the Eq transition of ^^Ge(n,n'e) could 
be simulated using Geant4.9.2. This makes the Geant4.9.2 simulation more 
reliable for dark matter searches utilizing germanium detectors. 



4-2 Verification of the Monte Carlo Simulation 



To obtain reliable results from our Monte Carlo, it was essential to verify 
the simulation. ^*^Co, with an original radioactivity of 1.0 /iCi, was used to 
validate the simulation. The ^°Co source was mounted on the center of the Ge 
detector cap as shown in Fig. [2] (right). The corresponding geometry used in 
the Monte Carlo simulation is shown in Fig. [7j Since the germanium detector 
is aged, the Monte Carlo geometry was modified to include a dead layer of 
thickness 2.5 mm (see Fig. ui left). We also implemented a smearing process 
by applying energy resolution to the spectrum of deposited energy in the 
active germanium. When fitting the energy resolution, these four peaks were 
identified: 1173 keV, 1332 keV, 2506 keV and 662 keV. The first three peaks 
are gamma rays from the ^"Co source and the last peak is a gamma ray from 



137 



Cs source. The best fit function is presented in Eq. 12 



0.2017 X J0.3+ (0.001E)(-i-96) + (0.001E)(-3-436)^ 



(12) 



where E is the energy in keV. Under the square root, the first term, 0.3, is 
the percentage of energy resolution, for a new Ge detector; the second and 
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Fig. 8. Percentage of relative energy resolution as a function of energy. 

third terms are due to noise and the age of the Ge detector. The percentage 
of relative energy resolution as a function of energy is plotted in Fig. |8] Using 
this fitted function, we applied the energy resolution to the range 500 keV ~ 
2510 keV and obtained agreement between the experimental data and Monte 
Carlo simulation as shown in Fig. |9| This validates that our Monte Carlo can 
be reliably used for gamma-ray simulation. 



5 Monte Carlo Simulation of Nuclear Recoils 



Using the modified Geant4.9.2 package, we performed simulations to verify 
the accuracy of the Lindhard and Barker-Mei models with an AmBe neutron 
source. The experimental setup for our Ge and Nal detectors, as shown in 
Fig. [2| was simulated according to the dimension and material information 
provided by the manufacturer pi]. Since the ^^Ge(n, ri'e) reaction causes a 
quasi-triangular shape in the data, we expect to see this same feature in our 
Monte Carlo simulation if the Lindhard and Barker-Mei models are accurate. 
In order to provide an accurate simulation, the geometry, as well as the AmBe 
source generator, need to be implemented correctly. 
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Fig. 9. Comparison between MC and data in energy deposition spectrum. 
5.1 AmBe Neutron Generator 



Because the AniBe source was placed very close to the germanium detector, 
gamma rays emitted from the source have a greater chance of entering the 
germanium detector and causing contamination in the region of interest (670 
keV ~ 760 keV). Thus, in the simulation it is necessary to account for all 
potential gamma rays emitted from the AmBe source. 



Two reactions, shown in Eqs. [13] and [14], occur in the AmBe source: 

2^Mm -^2^^ Np + a + -f (13) 

a+^5e^^^C + n + 7. (14) 



Eq. 13 shows the decay of ^^^Am to ^^"^Np, which causes the emission of alpha 



particles and gamma rays. Eq. [14] shows the reaction between an alpha particle 
and ^Be. In Eq. 14 , the energy of gamma rays emitted depends on the resulting 
state of ^^C, which is 4.443 MeV (for the 1st excited state), 7.65 MeV (for 
the 2nd excited state) and 9.64 MeV (for the third excited state). Using the 
gamma ray energy and recoil energy of ^^C, the resulting neutron energy can 



be calculated by applying energy conservation to Eq. 14 
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Fig. 10. Simulated neutron energy spectrum from AmBe neutron source. 



Fully absorbed gamma rays (from Eq. 14 ) will not be in the region of inter- 
est (670 keV ~ 760 keV) since their energies are at a few MeV. However, 
the Compton continuum of their interaction can contribute to the region of 
interest. In addition, gamma rays (from Eq. 13) with energies 26.34 keV (a 
branching ratio of 2.^ 



/o), 59.54 keV (a branching ratio of 35.9%) and 722.01 
keV (a branching ratio of 0.000196%) [35] can contribute to the region of in- 
terest by occurring in coincidence with the Eq transition (26.34 keV and 59.54 
keV) or by becoming fully absorbed (722.01 keV). We generated these three 
gamma rays in our Monte Carlo simulation along with gamma rays at ener- 
gies of 4.443 MeV, 7.65 MeV, and 9.64 MeV caused by the transitions in the 
excited ^^C nucleus. 



The simulated neutron energy spectrum from the AmBe neutron source is 



shown in Fig. 10 The spectrum agrees with the prediction from Marsh et. 
al. [36] . After this validation, the AmBe source generator was implemented in 
the Monte Carlo simulation. 



5.2 Verification of the Eq Transition in Simulation 



The simulated results are presented as a spectrum of energy deposited in the 
germanium detector after the application of smearing, which is the process of 
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Fig. 11. The spectrum of energy deposition in the germanium detector after smear- 
ing. 

accounting for energy resolution in a germanium detector. 

Gaussian and Moyal distributions with fitted parameters from our data (Sec- 
tion 3) have been incorporated into the smearing of different peaks in the 
region of interest. We used the model proposed by Lindhard et al. [3] (/c = 
0.159) to determine the ionization efficiency for a germanium detector. Fig. 11 
shows the simulated energy deposition spectrum after smearing. The £"0, 691.6 



keV, transition is visible, as shown in Fig. 11 Appearance of this peak indi- 



cates that the internal conversion process has been successfully implemented 
in our Geant4.9.2. 



6 Monte Carlo Shape Analysis 



After successfully simulating the AmBe neutrons in our Geant4.9.2, we col- 
lected the energy deposited in the germanium detector. This collected en- 
ergy spectrum was compared with our data after the application of smearing 



(Eq. 12). Two models utilizing ionization efficiency functions were applied to 
the Monte Carlo spectrum: Lindhard k = 0.159 and the Barker- Mei model. 
Normalization was applied to the energy range 670 to 760 keV. By overlaying 
the Monte Carlo energy spectrum and the collected data, we were able to per- 
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Fig. 12. Comparison of our Monte Carlo simulation utilizing the Lindhard (k = 
0.159) and Barker-Mei models and our collected data after background subtraction. 

Table 1 

The percent difference between the collected data and the Monte Carlo simulations 
with two models. E„r is the average visible nuclear recoil energy corresponding to 
each bin used in the analysis. 



Energy (keV) 


696 


702 


708 


714 


720 


726 


732 


738 


744 


Lindhard (k = 0.159) (%) 


1.4 


3.9 


3.1 


1.0 


0.8 


2.1 


2.8 


3.0 


2.9 


Barker-Mei (%) 


1.0 


3.2 


2.3 


0.6 


0.9 


1.8 


2.2 


2.1 


1.8 


E„, (keV) 


4.4 


10.4 


16.4 


22.4 


28.4 


34.4 


40.4 


46.4 


52.4 



form a shape analysis on the characteristic Eq transition. This is demonstrated 
in Fig. fT2l 



The shape analysis was performed on a bin-to-bin basis by comparing the data 
and two Monte Carlo simulations, which correspond to the two models. The 
difference in shape between the data and the two Monte Carlo simulations 
was calculated using (^ata+MC)/2 " Table 1 shows the results of this comparison 
in which the largest difference is shown to be less than 4%. 



Unfortunately, we were not able to include the neutron capture lines of even 
nuclei in our Geant4.9.2 package. Thus, at 708.2 keV and 747.7 keV, there is 
a discrepancy between the collected data and Monte Carlo simulations due to 
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the neutron capture of ^'^Ge(n, 7) [37]. There is also inconsistency in the energy 
range 650 - 680 keV that is hkely due to ^^^In(?2,7) and ^°^Pb(n, n'7) 
processes, which are not included in the Geant4.9.2 package. 



7 Conclusion 



We took measurements using an AmBe neutron source, incident on a germa- 
nium detector, for a total of 35 days. The characteristic quasi-triangular shape 
located at 691.6 keV represents the Eq transition of germanium-72. All peaks 
in the region of interest have been identified, thereby confirming source related 
and environmental backgrounds. The unique quasi-triangular shape induced 
by neutrons can be described using the Moyal distribution with the fitted 
parameters R = 0.06 and Empv = 696 keV, after background subtraction. Uti- 
lizing these parameters, we derived the most probable value for nuclear recoils 
as 4.4 keV electronic equivalent energy, which corresponds to a nuclear recoil 
energy of 17.5 keV. The average path length in the germanium detector for 
17.5 keV nuclear recoils was approximately 0.014 fim. 

A Monte Carlo simulation employing a corrected Geant4.9.2 package was mod- 
ified to duplicate the same experimental setup and AmBe neutron source. The 
Lindhard {k = 0.159) and Barker-Mei models were used to apply ionization 
efficiency to the energy spectrum and were compared to the experimental 
measurements. A bin-to-bin shape analysis was performed, and the difference 
between the measurements and two models were calculated. We obtained a 
percent difference that was less than 4% for the collected data and two Monte 
Carlo simulations. Because of this agreement, we calculated the most probable 
values for visible nuclear recoils (as shown in Table 1). Therefore, the Lindhard 
model (with k = 0.159) and Barker-Mei model can be used to determine the 
ionization efficiency in germanium detectors for 1 to 100 keV nuclear recoil 
energy. 
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